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We consider the two-band double exchange model for manganites with Jahn- Teller (JT) coupling 
and explore the suppression of the ferromagnetism because of the JT distortion. The localized spins 
of the t2g electrons are represented in terms of the Schwinger bosons, and two spin-singlet Fermion 
operators are introduced instead of the Cg electrons' operators. In terms of the new Fermi fields the 
on-site Hund's interaction is in a diagonal form and one accounts for it exactly. Integrating out the 
spin-singlet fermions, we derive an effective Heisenberg model for a vector which describes the local 
orientations of the total magnetization. The exchange constants are different for different space 
directions and depend on the density n of eg electrons and JT energy. At zero temperature, with 
increasing the density of the eg electrons the system undergoes phase transition from ferromagnetic 
phase (0 < n < nc) to A- type antiferromagnetic phase (ric < n). The critical value Uc decreases as 
JT energy is increased. At finite temperature we calculate the Curie temperature as a function of 
electron density for different JT energy. The results show that JT coupling strongly suppresses the 
spin fiuctuations and decreases the Curie temperature. 
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Jahn- Teller (JT) effect is related to systems with de- 
generated electronic states [H- The importance of the 
JT coupling for manganites was first discussed in 
with regard to the colossal magnetoresistance. The 
most widely studied representatives have chemical for- 
mula Rei-xA^MnOs, where Re is rare earth such as La 
or Nd, and A is a divalent alkali such as Ca or Sr. The 
important electrons in these compounds are Mn d elec- 
trons. They have five degenerate levels The crystal 
environment results in a particular splitting of the five 
d-orbitals {crystal field spliting) into two groups: the eg 
and <2g states. The electrons from the eg sector form 
a doublet, while the t2g electrons form a triplet. The 
population of the t2g electrons remains constant, and the 
Hund rule enforces alignment of the three <2g spins into a 
S = 3/2 state. Then, the t2g sector can be replaced by a 
localized spin at each manganese ion, reducing the com- 
plexity of the original five orbital model. The only im- 
portant interaction between the two sectors is the Hund 
coupling between localized t2g spins and mobile Cg elec- 
trons. The oxygens surrounding the manganese ion read- 
just their locations creating an asymmetry between the 
different directions. This effectively removes the degen- 
eracy of the Cg electrons' states. The lifting of the de- 
generacy due to the orbital-lattice interaction is called 
Jahn- Teller effect. 

The interaction between the electrons and phonons 
is unusually strong and leads to a wide range of 
striking physical phenomena. Changing the Cg elec- 
trons' concentration produces a variety of phases, which 
may be characterized by their magnetic, transport 
and charge-ordering properties The manganites 

Lai-xCaxMnOs have attracted interest due to their 
colossal magnetoresistance. The phase boundary be- 
tween ferromagnetism and paramagnetism, in these ma- 
terials, also separates a low temperature metallic phase 



from a high-temperature insulating phase. At temper- 
atures below Curie temperature T < Tc the resistivity 
is relatively low and increases as T is increased, whereas 
at T > Tc the resistivity is very high and (for most x) 
decreases as T is increased. The magnetoresistance for 
T K,Tc can be very large [3l ■ 

The double exchange model with JT coupling is a 
widely used model for manganites. The procedures fol- 
lowed to obtain the essential features of the model are 
different: numerical studies 0,01, Dynamical Mean-Field 
Theory (DMFT) [21,181,1^, ab initio density-functional cal- 
culations [13] , and analytical calculations 0, H, [O, ■ 
In spite of the common conclusion that JT coupling sup- 
presses the ferromagnetic state, the results are quite dif- 
ferent and do not match the experimental results. For 
example the calculated Curie temperatures are two and 
even three times larger then the experimentally mea- 
sured. Because of that it is important to formulate theo- 
retical criteria for adequacy of the method of calculation. 
In our opinion the calculations should be in accordance 
with the Mermin- Wagner theorem [l^. It claims that 
in two dimensions there is no spontaneous magnetization 
at nonzero temperature. Hence, the critical tempera- 
ture should be equal to zero. We employ a technique 
of calculation which captures the essentials of the 
magnon fluctuations in the theory, and for 2D systems 
one obtains zero Curie temperature, in accordance with 
Mermin- Wagner theorem. The physics of the ferromag- 
netic manganites near the Curie temperature is domi- 
nated by the magnon fluctuations and it is important to 
account for them in the best way. 

The present paper is focused on the influence of the 
JT distortion on the ferromagnetism of manganites. To 
model the manganites we employ the Hamiltonian H = 
Hde + Hei-ph- The first term describes the hopping of 
e„ electrons and the Hund interaction between the spin 
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Si of the Cg electron and the locaUzed t2g spin 

HdE = ^ tfi,C^i^Ct+sLl'a - 2Jif ^ • Si (1) 
iall'a i 

where c^^ and cua are creation and annihilation opera- 
tors for Cg electron with spin a in d^2_y2(d^^_^2) orbital 
at site i, and a is the vector connecting nearest -neighbor 
sites. For the cubic lattice, the hopping amplitudes be- 
tween I and V orbitals along the x, y, z directions are: 

tla = V3^:l, = ^/3^L = 34«b = ^ (2) 

The second term in Eq.([T|) is the Hund interaction be- 
tween the spin s.^ of the Cg electron and the localized <2g 
spin with sj" = 1/2 E cJ^^rtC,,^, where cr^,cry,cr"^ 

laf} 

are Pauli matrices, and the Hund's constant is positive 
{Jh > 0). 

The Hei-ph Hamiltonian models the coupling of Cg 
electrons to the lattice distortion 

k 

Hel-ph = 5^1 (<52'^" + Q3^rz^) + ^Y. (^2z + QI) (3) 

i i 

where T^i := ^ (c^aCzba + c^a'^^a^) and r^, = 
(c^aCiaa " c'l^^a'^iba) ■ In equation Q g is the electron- 

phonon coupling constant, while Q2i and Q^i are JT 
phonon modes. The second term in H^i-ph is the usual 
quadratic potential for distortions with constant k. The 
important energy scale of the phonon-electron interaction 
is the static JT energy Ejt = g^/(2fc). 

One can represent the spin operators Si of the local- 
ized t2g electrons in terms of Schwinger-bosons (lySia, (^J^) 
= I'^ta^ap'fiP^ V'taVia = 2s. By meaus of the 
Schwinger-bosons we introduce spin-singlet Fermi fields 

^^{r) = (r)c,,„(r) (4) 

V zs 

*»f W - ^ [mir)ca2{T) - f^2{T)cal{T)] (5) 

and write the spin of the Cg electron and the total spin of 
the system S'°* = S^ 4- in terms of the singlet fermions 
[3]. Further, we average the total spin of the system 
in the subspace of the singlet fermions A and B. The 
vector Mi = (S'°') / identifies the local orientation of 
the total magnetization. Because of the fact that ^29- 
electron spin is parallel with Cg-electron spin we obtain 

M, = f S, with M = S+^ Eii^ti^^fi - ^u^^i))f- 

I 

Now, if we use Holstcin-Primakoff representation for the 
vectors Mi (a+ , a) with M as an " effective spin" of the 
system (M^ = M^), the bose fields Oi and are the 
true magnons in the system. 



An important advantage of working with singlet 
fermions is the fact that in terms of these spin-singlet 
fields the spin-fermion interaction is in a diagonal form, 
the spin variables (magnons) are removed, and one ac- 
counts for it exactly. The theory is quadratic with respect 
to the spin-singlet fermions and one can integrate them 
out to obtain the free energy of fermions as a function 
of the magnons' fields a^, ai. We expand the free energy 
in powers of magnons' fields and keep only the first two 
terms. The first term -f/o, which does not depend on the 
magnons' fields, is a free energy of Fermions with spins 
of localized t2g electrons treated classically. We fix the 
model parameters and consider this term as a function 
of the JT distortion modes independent on the lattice 
sites . The numerical calculations shows that the func- 
tion depends only on vQl + Of ^^^^d set Q3 — 0. The 
physical value of the JT distortion is the value at which 
Ffo has a minimum. In this way we obtain the distortion 
as a function of the density of Cg electrons for different 
values of JT energy and fixed Hund's coupling. We fix 
the hopping parameter t = 1 to set the energy unit. The 
results for the renormalized distortion Q = gQ2 as a 
function of charge carrier density n are plotted in Fig. 1 , 
for different values of the JT energy Ejt and Jh = 15. 
The figure ([T]) shows that JT distortion appears at crit- 
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FIG. 1: (Color online) The renormalized distortion Q = gQ2 
as a function of density n for Jh ~ 15 and differnt values of 
the JT energy Ejt- Inset: The density n* , at which the JT 
distortion appears, as a function of JT energy Ejt- 

ical value of the charge carrier density n* and increases 
as density n is increased. The inset demonstrates that 
n* decreases and approaches zero as JT energy Ejt is 
increased. 

The second term in the Fermion free energy is 
quadratic with respect to the magnons' fields af,ai and 
defines the effective magnon Hamiltonian in Gaussian ap- 
proximation. 

i^cff = J2p'' ("-i^' + «i+aa»+a - a,^a»+a - a+^^a,) (6) 

ia 

In equation ^ are spin stiffness constants which de- 



FIG. 2: (Color online) Spin stiffness constants as a functions of density n for Jh = 15, Ejt = 1.73 (left), Ejt = 2 (middle) 
and Ejt = 4 (right). The vertical dash lines correspond to the density n* at which the JT distortion appears. 



pend on the space directions a. They are calculated 
at zero temperature, for fixed Hund's coupling, JT en- 
ergy, charge density, and JT distortion determined above. 
The calculations follow the technique developed in p^ . 
Based on the rotational symmetry, one can supplement 
the Hamiltonian ([6]) up to an effective Hcisenbcrg like 
Hamiltonian. written in terms of the vectors Mi 

i?cff = -^J"M, •M.+a (7) 
ia 

where — p'^/M. The ferromagnetic phase is stable 
if all effective exchange coupling constants are positive 
J*^ > {p'^ > 0). If one of them is negative, for ex- 
ample < {py < 0), and the others are positive 
> 0, > (p^ > 0,p^ > 0), the stable state is 
A- type antiferromagnetic phase which has planes {x, z) 
that are ferromagnetic (parallel moments), with antifer- 
romagnetic (antiparallel) moments between them. The 
spin-stiffness constant, as a function of charge carrier 
density, is depicted in Fig. ([2]) for Jh = 15, and three 
different values of JT energy, Ejt ~ 1.73, Ejt ~ 2 and 
Ejt = 4. The vertical dash lines correspond to the den- 
sity n* at which the JT distortion appears. The figure 
on the left illustrates in the best way the impact of the 
JT distortion on the spin stiffness constants. The ap- 
pearance of the distortion at n* is accompanied with a 
change of the slopes of the curves. The distortion splits 
the p*' (blue) and p^ (red) lines, and starts to de- 
crease. At critical density ric, p^ becomes equal to zero 
and the system undergoes a transition from ferromag- 
netic phase to A-type antiferromagnetic phase. The two 
other figures shows that spin stiffness constants decrease 
when JT energy increases and the critical density Uc de- 
creases too. As the spin stiffness constants are a measure 
for the magnon fluctuations in the ferromagnetic phase 
we conclude that JT distortion suppresses the magnon 
fluctuations. 

The most evident consequence of this suppression is the 
Curie temperature (Tc) decreasing. To calculate Tc we 
utilize the Schwinger-bosons mean- field theory [l^. We 
represent the vector Eq.([7]) by means of Schwinger 

bosons {cj)iaAta) 

= ^ E ^o.<P<I^^P = 2M (8) 

a/3 



Next we use the identity 

M, • M, = i (0+ (c^tp'f'^p) - I (0+ 0») (^^/s^j/s) 

(9) 

and rewrite the effective Hamiltonian in the form 

Hcff = E J'^i^U^+'^o.) {q^t+.Ap) (10) 

where the constant term is dropped. To ensure the con- 
straint we introduce a parameter (A) and add a new term 
to the effective Hamiltonian ([TU)) . 

HeS = i/eff + A E {<f>t'P^a - 2M) (11) 
i 

We treat the four-boson interaction within Hartree-Fock 
approximation. The Hartree-Fock hamiltonian which 
corresponds to the effective hamiltonian reads 

Hh-F = ^ E J'^^id+aUiS+a + ^E {'t'ta't'i'y ~ ^M) 
ia 

where Ui,i+a (wi.i-i-a) are Hartree-Fock parameters to be 
determined self-consistently. We are interested in real 
parameters which do not depend on the lattice sites, but 
depend on the space directions Mi,i+a = Ui i+a = "a- 
Then in momentum space representation, the Hamilo- 
nian p2|) has the form 

i?H-F = yE"a^""2AA/Ar + E£fc0fe0fc, (13) 

a fe 

where TV is the number of lattice sites and Sk is the dis- 
persion of the 0fc-boson (spinon). The free energy of the 
theory with Hamiltonian -ffn-F is 

F=^-Y^ulJ^-2XM+^-^Y.^n[l-e-*), (14) 

a k 

where T is the temperature. The equations for the pa- 
rameters Wa and A are: dF/du^ = dF/dX = 0. 

To solve the system of four equations it is more con- 
venient to introduce a new parameter (p,) instead of (A): 
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FIG. 3: (Color online) Tc as a function of eg electron density 
n for Jh = 15 and different values of the JT energy. 



A = ^ {ug^J'^ + ^Ua). In terms of the new parameter the 

a 

(/)fe-boson dispersion is = ^ [ua.J^ (1 — cos fca) + M^^a] 

a 

and the theory is well defined for positive constants 
Ma > and ^ > 0. For high enough temperatures fJ-{T) 
and Uii(T) are positive, and the excitation is gapped. De- 
creasing the temperature leads to decrease of fi{T). At 
temperature Tc it becomes equal to zero n(Tc) = 0, and 
long-range excitation emerges in the spectrum. There- 
fore this is the Curie temperature. We set n = and 
obtain a system of equations for the Curie temperature 
Tc and 

2 ^ — ^ cos k^' 

fc e ^ ■» - 1 



N Tgi-2: "aP''(l-COsfc,) 

A: e " - 1 

The results for the Curie temperature Tc as a function 



of Cg electrons density n are plotted in figure ([3]), for 
Jh = 15 and different values of the JT energy. The up- 
per (black) dash line is a reference line which corresponds 
to the case without JT distortion. The vertical dash dot 
lines, on the left, correspond to the density n*, while the 
vertical dot lines, on the right, correspond to the critical 
density Uc- The appearance of the JT distortion at n* 
leads to a spitting of the reference curve and the curve for 
a system with JT distortion. The density n* decreases 
when JT energy increases and the ferromagnetic phase 
is strongly suppressed because of the suppression of the 
magnon fluctuations, which in turn leads to the decreas- 
ing of the Curie temperature. 

To illustrate our results we present in a table the maxi- 
mal Curie temperatures {T™°'^[K]), for different JT ener- 
gies {EjT/t) and the corresponding Cg-electron densities 
(n). To do this we have utilized that t = O.'&eV 
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We have used a large value for Hund's constant to bet- 
ter demonstrate the impact of the JT distortion on the 
ferromagnetism. Decreasing of Jh suppresses the ferro- 
magnetic phase, decreases the Curie temperature, and 
reduces the impact of the JT distortion on the ferromag- 
netism. For example, for Jn/t = 15 and absence of JT 
distortion we have Tc = 1092iir, while for Jn/t = 5 we 
obtain Tc = 738X. For non-zero distortion, Ejx/t = 2, 
the Curie temperatures are Tc = SISK and Tc = 620iir 
respectively. These results show that the reduction of the 
Curie temperature due to JT distortion depends on the 
value of Jh- 
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